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ABSTRACT 



Context. The onset of the asymmetry of planetary nebulae (PNe) is expected to occur during the late Asymptotic Giant Branch (AGB) 
and early post-AGB phases of low- and intermediate-mass stars. Among all post-AGB objects, the most heavily obscured ones might 
have escaped the selection criteria of previous studies detecting extreme axysimmetric structures in young PNe. 
i-C ■ Aims. Since the most heavily obscured post-AGB sources can be expected to descend from the most massive PN progenitors, 

these should exhibit clear asymmetric morphologies. High-resolution observations of these sources should reveal marked bipolar 
morphologies, confirming the link between progenitor mass and nebular morphology. 

■ Methods. We have obtained VlSlR-VLT mid-lR images of a sample of four heavily obscured post-AGB objects barely resolved 
^ ■ in previous Spitzer IRAC observations in order to analyze their morphology and physical conditions across the mid-lR. The images 
^ I obtained in four different mid-lR filters have been deconvolved, flux calibrated, and used to construct RGB composite pictures as well 

as color (temperature) and optical depth maps that allow us to study the morphology and physical properties of the extended emission 
of the sources in our sample. 

Results. We have detected extended emission from the four objects in our sample and resolved it into several structural components 
^ ' that are greatly enhanced in the temperature and optical depth maps. The morphologies of the sample, as well as their physical 

I conditions, reveal the presence of asymmetry in three young PNe (IRAS 15534-5422, IRAS 17009-4154, and IRAS 18454+0001), 

1/^ ■ where the asymmetries can be associated with dusty torii and slightly bipolar outflows. The fourth source (IRAS 18229-1127), a 

■ possible post-AGB star, is better described as a rhomboidal detached shell. 
Conclusions. The heavily obscured sources in our sample do not show extreme axisymmetric morphologies. This is at odds with the 
expectation of highly asymmetrical morphologies in post-AGB sources descending from massive PN progenitors, which is otherwise 
supported by observations of bright mid-IR unobscured sources. The sources presented in this paper may be sampling critical early 
phases in the evolution of massive PN progenitors, before extreme asymmetries develop. 



Key words, techniques: high angular resolution - Stars; post AGB -circumstellar matter - planetary nebulae: general 



X 



1 . Introduction ing mechanism of complex PNe i s a matter of intense de bate and 

H"^,., .. ,. is not vet completely understood (iBalick & Fran5l2002h 
^ . Evolved low- and mtermediate-mass stars (O.8M0 < M < 8M0) 

undergo heavy mass-loss episodes during the red giant and has been claimed that the short transition between the 

Asymptotic Giant Branch (AGB) phases. At the tip of the AGB, ^GB and post-AGB phase contams the clues to understand 

these stars will eject most of their stellar envelopes in a short *e transformation of the spherical_stellar envelope mto an 

time until they evolve into post-AGB stars, the immediate pre- asymmetric PN (Sahai & Traugei] [1998D. AGB stars are sur- 

cursors of planetai-y nebulae (PNe). PNe display an impressive rounded by thick and compact circumstellar envelopes rich in 

variety of morphologies, from spherical to the most complex dust; this makes their detection almost impossible at optical 

shapes; bipolai-s and multipolars with point-symmeti'ic struc- wavelengths, whereas their infrared emission is strong. As they 

tures and coUimated jets moving at high velocities. It has been evolve into the post-AGB phase, the temperature of the stel- 

suggested that asymmetric PNe are the rule, rather than the ex- l^r cores increases and the envelopes become optically thin, 

ception (Manchado et al. 2000; Schwarz et al. 1993). The shap- although it must be noted that there are sources at this stage 

and even young PNe that are significantly obscured because 

* Based on observations collected at the European Organisation for they preserve optically thick envelopes (e.g., IRAS 17347-3139, 

Astronomical Research in the Southern Hemisphere, Chile. Program: de Gregorio-Monsalvo et al.1 2004). Observational studies of the 

087.D-0367(A). reflection, thermal dust and ionized emission around proto- 
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Fig. 1. Near- and mid-IR spectral energy distributions (SEDs) of the sources in our sample. See lRamos-Larios et al. I (120091 120TI 
for further details on the different datasets used to built these SEDs. 



PNe and young PNe ISahai & Traugedl 19981: ISahai et alJl2007t 
lUeta et al.l 120001: iLagadec et al 1 1201 U ISahai et all 1201 Ih typi- 
cally detect the existence of extreme axisymmetric morpholo- 
gies, with highly collimated bipolar lobes and equatorial rings 
or torii. On the other hand, spherical and elliptical morphologies 
seem to be rare among the objects in this transition phase. 

Infrared surveys (e.g., 2MASS, IRAS, wdAKARI) have con- 
tributed to the study of objects in the late AGB and early post- 
AGB stages, providing the first insight s on this short transi- 
tion phase (e.g ., Garcfa -Lario et al.lll997 : Jimenez-Esteban et al. 
2006[ ICox et al.i ,201 ll) . Recently. iRamos-Larios et al., (,2009>. 
2OI2I hereafter RL09 and RL12, respectively) have inves- 
tigated a sample of heavily obscured post-AGB and PNe 
candidates selected according to their IRAS col ors and the 
lack of an optical counterpart (ISuarez et al ] l2006h using near- 
IR JHK and Spitzer IRAC GLIMPSE images, and MSX, 
AKARI, and IRAS photometric data. RL09 and RL12 have 
gathered a small sample of four objects that are resolved in 
Spitzer IRAC GLIMPSE images, namely IRAS 15534-5422, 
IRAS 17009-4154, IRAS 18229-1127, and IRAS 18454H-0001. 
The spatial resolution of these images (~2"), however, is not 
adequate to investigate the angularly small structures t ypical 
of sources in these evolutionary phases dSuarez et al.i 1201 it 
iLagadec et al ] l20Ilh . Mid-IR observations using the new gen- 
eration of ground-based telescopes provide a unique opportunity 
to resolve the extended emission detected in these four highly 
obscured post-AGB objects in their transition to the PN phase. 

In this paper we present mid-IR VISIR-VLT high angular 
resolution observations of these four sources. The images have 



been used to describe their morphologies, to analyze their in- 
nermost structures, and to derive color (temperature) and optical 
depth maps. We next describe the sample in Sect. 2, the obser- 
vations and data reduction in Sect. 3, and present the results in 
Sect. 4. The results are discussed in Sect. 5 and a short conclu- 
sion is provided in Sect. 6. 

2. The sample 

RL09 and RL12 reported the detection of extended emis- 
sion in Spitzer IRAC images of four post-AGB source 
candidates, namely IRAS 15534-5422, IRAS 17009-4154, 
IRAS 18229-1127, and IRAS 18454+0001. According to the 
spectral energy distribution (SEP) classification schem e of post- 
AGB sources introduced bv lvan der Veen et~an (Il989l) . the SEDs 
of these four sources shown in Figure[T]can be assigned to Type 
II for IRAS 15534-5422 (peak at -25 jum and gradual fall-off to 
shorter wavelengths, although we note that its SED also suggests 
a near-IR excess at 2-5 pm), to Type III for IRAS 18229-1 127 
and IRAS 18454+0001 (peak at -25 p\n and steep fall-off to 
shorter wavelengths with a plateau between 1 and 4 pm), and 
to Type IV for IRAS 17009-4154 (with a main peak at 25 pm 
and a secondary peak blue-wards). These SEDs are suggestive 
of two dust components, cold dust in the thermal IR and hot dust 
obscuring the central star in the near-IR. 

Besides the information provided by RL09 and RL12, 
there is no detailed study available in the literature for these 
sources. Based on their IRA S colors, IRAS 15534- 5422 is clas- 
sified as a PN candidate (IPreite-Martined ri988h . as recently 
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confirmed by iParker et~an (1201 2h by means of spectroscopy, 
whereas IRAS 18454+0001 is classified as a post-AGB star 
(Garci'a-Lario et al. 199 3)- We note th at the IRAS selection cri- 
teria for these sources dSuarez et alJ 1200 6) may overlap with 
those of Young Stellar Objects (YSOs), although RL09 noted 
that these four objects are not located near star forming re- 
gions, neither they have been classified as YSOs in the litera- 
ture. The lack of CO Hne emission in IRAS 15534-5422 and 
IRAS 17009-4154, and the detection of narrow (^^0.8 km s'^) 
CO l ine emission toward IRAS 18454+0001 dUrquhart et alJ 
l2008l RMS Survey) are inconsistent with a YSO nature. Based 
on the IRAS colors of the sample, we can also discard a possible 
symbiotic star nature, as these have typica lly values of the IRAS 
[I2]-[15] color -0.8 dKenvon et al.lll988l) . whereas the sources 
in our sample present values >2. Moreover, our sources do not 
present variability in near-IR observations (RL09 and RL12), 
thus we can disregard an eruptive nature, as for instance in cir- 
cumstellar shells a round luminous bl ue variable sources (e.g., 
IRAS 18576+034I. [Buemi et aLllIoToh . 

The detection of radio continuum emission in 
IRAS 15534-5422, IR AS 17009-4154, and IRAS 18454+0001 
dUrquhart etalJ l2008l) confirms that ionization is already 
present, implying that these three sources may have already 
entered the PN phase. Indeed, Bry emission from ionized 
material has been detected in IRAS 15534-5422 ( RL12) 
and in IRAS 17009-4154 dVan de Steene et al.1 l200tf ). For 
IRAS 15534-5422, this emission is found to be extended and to 
display a bipolar morphology (RL12). We are thus confident on 
the PN nature of IRAS 15534-5422, IRAS 17009-4154, and 
IRAS 18454+0001. 

As for IRAS 18229-1 127, the absence of filamentary diffuse 
emission in the mid-IR, which is characteristic of YSOs, and the 
similarities between its SED and that of IRAS 18454+0001 can 
be used as additional arguments for a post-AGB classification. 



3. Observations and data reduction 

High angular resolution observations (Program ID: 087. D- 
0367(A), PI : M.A. Guerrero) w ere obtained with the mid-IR im- 
ager VISIR dLagage et al.l2004b attached to the Cassegrain focus 
of Melipal (UT3) at the VLT. The sources have been observed 
through four different filters, PAHl (Ac-8.54 yum, AA-0A2 yum), 
SiC (Ac=ll-85 fim, A^=2.34 fim) and [Neii] (^1^=12.82 fim, 
A/l=0.21 ^m) in the N band, and the Ql filter (/lc=17.65 yum, 
A/l=0.83 fim) in the Q band. The observations of standard stars 
were performed after every science ObsBlock to coiTect for PSF 
artifacts and flux calibration purposes. The date of observation, 
the integration time, and the FWHM of the standard star ob- 
served subsequently are summarised in Table [1] The pixel scale 
is 0'.'075 and the field of view (FoV) is 19'.'2xl9'.'2. 

The data were taken using several exposures with short DIT 
(Detector Integration Time) depending on the flux of the source 
and the required S/N. The observation mode used was the so- 
called NORMAL with a perpendicular chop-throw of 8". The 
chopping and nodding standard technique was used to help in 
the removal of the background signal. In the NORMAL mode 
used for these observations, all frames taken at a chopping posi- 
tion are added immediately at the end of the exposure, resulting 
in a data cube of reduced size. The data reduction was carried out 
following standard procedures of Gasgano-VISIR pipeline (ver- 
sion 3.4.4), in which flat fielding correction, bad pixel removal, 
source alignment, and co-addition of frames are executed to pro- 
duce a combined image for each filter The resulting FoV is <8". 



In order to improve the spatial resolution of the raw im- 
ages, we deconvolved each of them using its PSF obser- 
vation according to two different deconvolution algorithms: 
Maximum likelihood (number of iterations >5) and Richardson- 
Lucy (number of iterations> 10). The precise algorithm and 
number of iterations was determined by the quality of the 
raw image, so that the spatial resolution of the decon- 
volved image was improved but no artifacts were intro- 
duced. The maximum likelihood algorithm was used to de- 
convolve the images of IRAS 17009-4154, IRAS 18229-1127 
and IRAS 18454+0001, whereas the Richardson-Lucy algorithm 
was used for IRAS 15534-5422. The final images are presented 
in Figure 121 

The final ima ges have been fl ux calibrated using fluxes of 
the PSF stars (iCohen et alJll999\ For the flux calibration we 
have performed aperture photometry (Table |2]). We note that the 
angular extent of IRAS 15534-5422 and IRAS 17009-4154 is 
similar to the final FoV of 8". As the background aperture may 
include emission from these objects, their fluxes listed in Table|2] 
should be regarded as lower limits of the real fluxes. 

4. Results 

The VISIR- VLT images presented in Figure |2] reveal with 
unprecedented detail the extended emission from these four 
sources. These images also disclose notable variations among 
the morphological features shown in images obtained through 
different filters which we interpret to be related to the spectral 
features registered by these filters. The emission in the PAHl 
filter may include the PAHl feature at 8.6 jum associated with 
C-rich dust, as well as thermal continuum emission. The broad 
SiC filter includes thermal dust emission, silicate features, the 
PAH2 feature at 11.3 //m, and emission lines such as [Ne ii] and 
[S iv]. The [Neil] filter is fine tuned to the [Neii] emission from 
ionized gas, but it may also include the contribution from dust 
emission. Finally, the Ql filter maps the continuum emission of 
the thermal dust at 17.7 jum. 

Considering that the thermal emission of the dust represents 
the major contribution to the emis sion in the mid-IR we have 
applied the procedure described bv lDaval et alj dI998h . we have 
generated color maps (also known as temperature maps) using 
pairs of flux calibrated images in different wavelengths for the 
four objects in our sample (Figure [3]). This procedure relies on 
the relation between the intensity of the thermal dust continuum 
at a given wavelength (/^) with the temperature and optical depth 
under the assumption of the optically thin emission of the warm 
dust component. As the extinction produced by cool dust com- 
ponent along the line of sight can be considered to be similar 
at the mid-IR wavelengths of the different filters used in this 
work, we can further assume that its effects on the ratio maps 
are negligible. Since variations in the temperature produced by 
the emission properties of the dust appear as color variations in 
the maps, the temperature of the dust can be approximated by 
the following expression: 

^ 1.44x10^(1/^2- IMi) ^^ 
ln[(hJIM/A2)^] 

where Ai and A2 are the two wavelengths used for the estimation 
of the temperature, and and are the measured intensities at 
these two wavelengths. This method assumes that the variation 
of the emissivity (0 can be represented by a power law (Q ~ 
A'") and it has been proven to provide a good estimate of the 
spatial distribution of the cold and warm dust across the source 
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Fig. 2. Deconvolved VISIR images of the objects in our sample in the different filters observed in the N and Q bands. North is 
up and east to the left. The images of IRAS 18454+0001 are overlaid with contours at 30, 40, 50, 60, 70, 80, 90 and 100% of 
the peak intensity to highlight its morphological features. The frames at the bottom are color-composite RGB pictures using the 
filters with the colors of the corresponding labels. The bar at the bottom of each RGB image and at the bottom of the Ql image 
of IRAS 18454+0001 represents 1". The scale of the color-composite picture of IRAS 18454+0001 is different from that of the 
individual images to show in the former the faint halo detected in the Ql band. 
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Table 1. Observations log 



Object 


or S 


Date of 


Integration time 


Standard star 




FWHM 






(J2000.0) 


observation 


N band 


Q band 




PAHl 


SiC 


[Nen] 


Qi 








(s) 


(s) 




(") 


(") 


(") 


(") 


IRAS 15534-5422 


15''57'"21'';il -54°30'46'.'4 


2011-05-25 


600 


900 


HD 133550 


0.28 


0.37 


0.35 


0.46 


IRAS 17009-4154 


17hQ4m29;5 -41°58'38'.'9 


2010-05-24 


300 


450 


HD 163376 


0.31 


0.36 


0.38 


0.48 


IRAS 18229-1127 


18''25"'45^.0 -ir25'56" 


2010-06-14 


600 


600 


HD 169916 


0.34 


0.51 


0.42 


0.49 


IRAS 18454+0001 


18H8'"0P.5 +00°04'47" 


2010-06-22 


840 


1320 


HD 168723 


0.28 


0.35 


0.35 


0.49 



Table 2. Estimated fluxes 



Object 


PAHl (8.54 ^m) 


SiC(11.9yL(m) 


[Neil] (12.8 /im) 


Ql (17.7 yum) 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


IRAS 15534-5422 


181.6±0.1 


768.5+0.1 


1461.9+0.1 


1891.9+0.2 


IRAS 17009-4154 


9054.5±0.5 


4212.2+0.1 


3280.2+0.2 


50674+5 


IRAS 18229-1127 


659.78+0.02 


270.25+0.04 


270.79+0.03 


4101.8+0.1 


IRAS 18454+0001 


184.34+0.03 


167.61+0.02 


225.17+0.03 


3096+1 



studied iD aval et al.'l998V Meixner et all 1 9991: lUeta et al.l200U 
iLagadec et al..2005i) . being helpful to reveal or enhance morpho- 
logical features. 

In order to derive the color maps shown in Figure[3] we have 
used the Ql image in conjunction with one N-band image, as the 
longer wavelength range provides a better leverage of the contin- 
uum slope which is determined by the thermal dust temperature. 
We have also avoided images in those filters that may present 
important contributions of spectral features others than the con- 
tinuum. Therefore, the SiC images have not been used, as they 
may include multiple contributions from emission lines and dust 
features. 

Once the value of the temperature is estimated, we can con- 
struct optical depth maps using the Planck function under the 
same assumption. Therefore, the optical depth at a certain wave- 
length A can be derived using the following expression: 



-In 



1 



Ba(T) 



(2) 



where I a is the intensity at the wavelength analyzed. Using these 
maps, we are able to estimate the column density variations 
across the nebula to assess whether the emission is optically thin 
at the observed wavelength range. 

We next describe in detail the morphology, temperature and 
optical depth properties of the individual sources. 

4.1. IRAS 15534- 5422 

IRAS 15534-5422 is resolved in at least two structures that 
are distinguished by their morphology and emission properties 
(Figure|2]i. A first structure is an elongated, bar-like feature with 
a size of ~3" oriented at PA 55°. The emission of this bar is de- 
tected in PAH, SiC, and Ql, and is much weaker in [Neii]. A 
second structure is an arc-like feature, ~5"x 3" in size, which 
is oriented perpendicular to the bar and is particularly prominent 
in [Neil], but weak in SiC and Ql and absent in PAHl. This 
arc-like extended structure is not closed and presents a rectan- 
gular shape at low intensity levels. The color composite picture 
of IRAS 15534-5422 (Figure |2]i suggests that the arc-like fea- 
ture traces ionized material, while the bar is mostly dominated 
by thermal dust emission. 

The color map of this source (Figure [3]), obtained using the 
[Ne ii] and the Ql images, reveals a range of temperatures across 
the nebula from 80+1 K to 270±1 K, with a mean value of 



~190 K. The color map greatly enhances the arc-like feature, 
which is the hottest component, whereas the bar, with lower tem- 
peratures, almost disappears in this map. The optical depth maps 
at different wavelengths (Figure HI, computed for a mean tem- 
perature value of 190 K, imply that the south-west tip of the bar 
represents the highest density zone in this object. The [Ne ii] op- 
tical depth map reveals a noticeable enhancement of the column 
density at the north-east tip of the bar, where it goes across the 
arc-like feature. 



4.2. IRAS 17009-41 54 

IRAS 17009-4154 has an elongated morphology toward the 
north-west with a size of 6'.'3x5'.'5 in size. There is an innermost 
equatorial enhancement that can be appreciated as an elongated 
ring-like structure (~3"x2") at PA~64° surrounding a central 
star which is detected in all VISIR bands (Figure|2]i. Two promi- 
nent arc-like features trace the edge of the outermost regions 
and are located at PAs 64°-244° with brightness increasing with 
wavelength. 

The color map of this source (Figure O shows values of 
temperature in the range from 150±5 K up to 390+5 K with a 
mean value of 210 K. The maxima in temperature are related 
with regions at the tips of the major axis of symmetry. The in- 
nermost ring-like structure suiTounding the central star seems to 
have two different temperature components: the south-west is 
the coolest region (~150 K), whereas there is a hot component 
located toward north-east with a temperature of ~350 K. The 
arcs are warmer (~250 K). The mean value of the temperature 
(r=210 K) has been used to derive the optical depth maps of 
IRAS 17009-4154 (Figure|4|. The regions with the highest col- 
umn densities are the arcs, followed by the ring-like structure, as 
shown in all optical depth maps. 



4.3. IRAS 18229- 11 27 

IRAS 18229-1127 has a rhomboidal clumpy envelope with an 
angular size of 2'.'5. In the PAHl, SiC, and [Neii] N-band im- 
ages (Figure |2]i, this envelope is dominated by two bright knots 
toward the north and north-west, and two weaker knots are no- 
ticeable toward the south and south-east. The south-east knot is 
not detected in [Ne ii] . Interestingly, these knots are singularly 
not apparent in the Ql image that shows a clearly hollow enve- 
lope, with a elongated cavity oriented along the north-south di- 
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RAS 18229-1127 - 8.6 / 17.7 microns IRAS 18454 + 0001 - 12.8 / 17,7 microns 




Fig.3. Temperature (color) maps derived using the equation 1 of IRAS 15534-5422, IRAS 17009-4154, IRAS 18229-1127 and 
IRAS 18454-18454. We note that the red patches at the borders of the color the maps of IRAS 17009-4154 and IRAS 18454+0001 
are artifacts caused by the reduced S/N ratio in the outermost regions of these sources. The central stars of IRAS 17009-4154 and 
IRAS 18229-1 127 have been masked to construct the color maps 



rection. The central star is detected, but only in the bluest PAHl 
image at 8.6 /im. The color-composite RGB picture of this object 
at the bottom of Figure |2] shows a reddened, hollow rhomboidal 
envelope with a spot brighter at shorter wavelengths at each cor- 
ner. 

We have derived the color map of this source using its R\H1 
and Ql images and we estimate a mean temperature of 170 K 
(FigurelS, with a range of 100 K< T < 370+1 K. 

The bright knots in the N-band filters are the warmest regions 
of this source, with T ^280 K, whereas the rest of the rhom- 
boidal envelope has the lowest temperatures (~200 K). The mean 
value of 170 K has been used to derive the optical depth maps of 
IRAS 18229-1127 (Figure H)). The knots represent the regions 
with the highest densities in the N-band, with the north one be- 
ing the most dense. These density variations are diminished in 
the Ql optical depth map, where the rhomboidal envelope has 
a more homogeneous density distribution. It is remarkable the 
difference between the optical depth maps of the N-band and the 
one of Ql: while the knots are the densest zones from 8.6 to 12.8 
/im, the peak of the density at 17.7 //m is toward east, where not 
a single knot is present. 



4.4. IRAS 18454+0001 

This source has the smallest angular size among the objects in 
our sample. The source can be described as a round, r.'5 in size, 
clumpy envelope in the SiC, [Ne ii], and Ql images, but its mor- 
phology in the R4H1 image is singularly different (Figure|2]l. In 
this band, the image is dominated by a pair of bright knots along 
PA 121° that seem spatially coincident with brightness peaks in 
the Ql image, whereas the peaks detected in the SiC and [Neii] 
images are located at different positions. All images hint at the 
presence of a central cavity, but its orientation and extent varies 
with wavelength. At large spatial scale, the Ql image reveals 
an extended halo around the main nebula. The halo with a size 
of 5"x4" extends toward the north-east, as shown in red in the 
color-composite RGB picture at the bottom Figure |2] 

The [Neil] and Ql images of IRAS 18454+0001 have been 
used to estimate the mean value of temperature of 80 K in a range 
from 50±1 K to 120+1 K (Figure|3]l. The temperature variations 
are notably flat over the nebula, with the hottest temperature 
(>100 K) at the pair of bright spots seen in [Neii] and regions 
of high temperature (^90 K) at the bright PAHl and Ql knots. 
The emission of the envelope seems to be traced by a cooler 
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dust component. Interestingly, the optical depth maps (Figure|4]l 
show that the highest density regions are associated with the 
knots seen in PAHl and in Ql, while the knots associated with 
the [Ne ii] emission peaks represent zones of lower column den- 
sity. There is a remarkable difference in the value and location 
of the peaks in the the optical depth maps of IRAS 18454+0001 
in PAHl and Ql with respect to those in SiC and [Ne ii], but 
we note this may be an artifact produced by the low value of the 
temperature estimated for this source (7=80 K) that has been 
used to built its optical depth maps. 

5. Discussion 

5.1. Physical structure of the sources 

We have obtained mid-IR high spatial resolution images of three 
PNe and one post-AGB source candidate using VISIR at the 
VLT. These images reveal notable morphological differences 
among the sources in our sample. The shell-like morphology 
and detection of a central star in three of these sources are not 
expec ted for the diffuse mid-IR emission from compact H ii re- 
gions (iTakahashi et al.l2 000). Only IRAS 15534-5422 may look 
like a compact H n regi on, but its PN natur e is well established 
by optical spectroscopy dParker et al.ll20I2|) . 

Temperature and optical depth maps built from these images 
are proven to be powerful tools to enhance morphological struc- 
tures and to reveal new ones in these objects. In the following we 
will discuss the main characteristics of the individual sources as 
well as the characteristics that are shared by them all. 

IRAS 15534-5422 is a young PN with two distinct structural 
components. The elongated innermost structure is dominated by 
thermal dust continuum and it represents the coldest and dens- 
est regions in this PN. This structure is surrounded by a warmer, 
ionized extended envelope of lower density. Whereas the color- 
composite picture of IRAS 15534-5422 in Figure |2] seems to 
imply that the elongated inner structure is a bipolar nebula sur- 
rounded by the ring-like envelope, the physical properties of the 
different components favor an alternative, more solid interpre- 
tation. The innermost structure is rather a high-density, dusty 
"torus", whereas the apparent envelope can be associated to a 
pair of ionized bipolar lobes perpendicular to that "torus". This 
interpretation is consistent with th e narro w -band near-IR images 
presented by Ramos-Larios et alj (l20I2h and with the near-IR 
excess noticeable in its SED (Fig.[l)- 

IRAS 17009-4154 seems to have an slightly bipolar or el- 
liptical morphology. The bright arcs can be interpreted as ion- 
ized bipolar extensions of an elliptical shell sculpted by a SiC 
and Ql bright, dusty equatorial ring. The dust ring is the coldest 
structure in this source, whereas the bright arcs are warmer The 
temperature is the highest at the northwest and southeast tips of 
the lobes. 

Based on the properties of the SED of IRAS 18229-1127 
and on its similarity with that of the PN IRAS 18454-1-0001 
(Fig. [TJ, we favor its classification as a post-AGB source and 
discard a YSO nature. This source shows a peculiar clumpy and 
dusty rhomboidal envelope. This rhomboidal appearance is en- 
hanced by the knots located toward the northwest, the densest 
region of this source, and southeast knots. The rhomboidal mor- 
phology of I RAS 18229-1127 resembles the shape of t he PN 
BD+303639 (iLaeadec et alllloTltlAkras & Steffenll20I2l) . 

IRAS 18454-1-0001 is a compact PN displaying both a spher- 
ical AGB shell and a dusty torus, most noticeable in the PAHl 
(8.6 jum) and Ql (17.7 fim) images. The temperature map sug- 
gests the presence of relatively hot regions along the direction 



perpendicular to this torus. Alternatively, these regions may be 
dominated by ionized gas unveiled by its [Ne ii] emission. The 
torus shows density enhancements in the dust-dominated PAHl 
and Ql filters corresponding with the high density knots along 
the south-east and north-west direction. The presence of a high 
density torus, but a spherical AGB envelope, is intriguing. We 
might witness in this source the early phase on the shaping of a 
bipolar PN, before the spherical geometry of the AGB wind is 
disrupted by a fast stellar wind collimated by a toroidal density 
enhancement as proposed by the G eneralized Interacti ng Stellar 
Wind model of PN shaping (GISW lBalick et al.|[T987l) . 

The morphological features detected in direct images are no- 
tably enhanced in the temperature and optical depth maps, as 
these distinguish between structures of same brightness, but dif- 
ferent physical conditions. Furthermore, the physical character- 
istics of the structures detected for each source are generally con- 
sistent. The innermost dusty rings or torii of IRAS 15534-5422, 
IRAS 17009-4154, and IRAS 18454+0001 have lower temper- 
atures (T -100 K) than the elongated hotter (T >150 K) bipolar 
regions. 

The optical depth maps of IRAS 15534-5422, 
IRAS 17009-4154, and IRAS 18229-1 127 peak at similar 
values for each of these sources. This implies that, there are 
no variations in the density derived from the different filters 
for each of these sources, and thus that the mid-IR emission of 
these source is optically thin. 

5.2. Comparison with the morphology and physical structure 
of other evolved objects 

Our mid-IR VISIR- VLT images have unveiled a wealth of 
structural components in the four sources in our sample. In 
short, these sources can be described as a detached rhom- 
boidal shell (IRAS 18229-1127), a spherical shell with strong 
evidences of a dusty torus orthogonal to an ionized bipo- 
lar flow. (IRAS 18454+0001), and two mild bipolar or el- 
liptical sources (IRAS 15534-5422 and IRAS 17009-4154). 
Contrary to other st udies of sources in the transition towards the 
PN phase (Sahai & Trauger 19981: lUeta et al.ll2000l: ISahai et"al] 
2007, 201 1; Lagadec et al. .201 1) , we do not find evidence of 
extreme axisymmetric bipolar or multipolar morphologies. In 
this respect, we note that we have not imposed to the sources 
in our sample any of the restrictions of previou s studies, that re- 
quired the source s to be detected in the optical ( Ueta et aDl2000l : 
I Sahai et a ni20 1 1[). or to have fluxes at 12 jum greater than 10 Jy 
(Lagadec et al] |20111) . These selection criteria may have intro- 
duced notorious biases in the morphological output of the sam- 
ple. 

In this sense, the lack of restrictive selection criteria in 
our sample may probe a different population of sources in 
the transition between the late AGB and early post-AGB 
phases. The small number of objects in our sample does 
not allow us to drawn firm conclusions, however, the ob- 
served morphologies may be analyzed within the evolution- 
ary context of PN formation. At least three of the sources 
in our sample (IRAS 15534-5422, IRAS 17009-4154, and 
IRAS 18454+0001) have already reached the PN stage, although 
they are still mostly obscured at optical wavelengths. The signif- 
icant obscuration can be attributed to large amounts of circum- 
stellar material that has been previously ejected by their central 
stars, suggesting massive progenitors. Interestingly, the optical 
depth values computed for IRAS 18454+0001 are significantly 
larger than those of the remaining sources, whereas its tempera- 
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Fig. 4. Optical depth maps of the sources obtained using equation 2 for the four different wavelengths observed. Contours are in 
intervals of 10% of the intensity peak for IRAS 17009-4154, IRAS 18229-1127 and IRAS 18454+0001. For IRAS 15534-5422 
the contours are in intervals at 10% for [Ne ii] and SiC, intervals at 20% for Ql, and for the PAHl map the contours start at 70% of 
the intensity in intervals of 2.5% of the intensity peak. 



ture is lower. The thicker envelope of this source may be indica- 
tive of an even more massive progenitor. 

If the sources in our sample indeed descended from mas- 
sive AGB stars, we do not see the expected correlation 
betw een extreme bipolar morphology and rnassive progeni- 
tors ('Corradi&Schwa rzlll995t lUeta et al.ll2000t ISiodmiak et aP 
[2OO8). It can be argued that the sources in our sample are at a 
very early evolutionary phase, when asymmetries have not fully 
developed yet. For example, IRAS 15534-5422, with its dense, 
dusty torus and bipolar lobes, might be a bipolar PN in the mak- 
ing, whereas the high density torus of IRAS 18454+0001, en- 
closed within a spherical shell, can provide the seed for a bipolar 
PN. 

The use of color maps in the mid-IR domain has revealed 
a torus in IRAS 07134+1005 (Dayal et al. 1998) and hinted the 
prese nce of jets in Roberts 22 and in V Hya dLagadec et alj 
I2005h . The mid-IR size of these sources (~4") are similar to that 



of the sources in our sample, but the variations of the temperature 
across the sources studied in this work have shown structures 
that define the morphology more clearly. Compared with these 
previous analyses of color and optical depth maps, our study re- 
solves much more morphological details, specially if we note 
that the sources in our sample do not display extreme axisym- 
metric nebulosities. The high-spatial resolution achieved by our 
data is certainly required to resolve the extended emission of the 
obscured and small-sized or compact sources transiting this evo- 
lutionary phase. 



6. Conclusions 

We have observed four heavily obscured post-AGB sources 
and PN candidates with VISIR-VLT in three different N bands 
(PAHI, SiC and [Ne 11]), and one Q band (Ql). Three sources 
in our sample (IRAS 15534-5422, IRAS 17009-4154, and 
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IRAS 18454+0001) can be classified as young PNe, whereas 
only a preliminary classification as a post-AGB source is pos- 
sible for IRAS 18229-1127. 

The high-spatial resolution VISIR images have been used 
to investigate the extended emission and to study the spatial 
variations of the physical conditions (temperature and optical 
depth) of these sources. We are providing evidence of asym- 
metry in three young PNe: an innermost dust torus or ring em- 
bedded within an ionized spherical shell (IRAS 18454+0001), 
and two mild bipolar or elliptical sources with dusty rings 
(IRAS 15534-5422 and IRAS 17009-4154). 

Compared to previous works, the use of color and optical 
depth maps have proven much useful to reveal fine structural de- 
tails in a small sample of heavily obscured sources reaching the 
PN phase, confirming the usefulness of mid-IR high resolution 
observations for the study of this short evolutionary phase. 

Our imaging study confirms that asymmetry is present in 
heavily obscured post-AGB stars and young PNe but, contrary 
to previous studies focused in these evolutionary phases, no ex- 
treme axisymmetric morphologies are found. These previous 
studies may be biased towards the mid-IR and optical brightest 
sources, missing critical early phases of the post-AGB evolution 
of the most massive progenitors. Our study, although based on 
a small sample, may be yielding important clues on the onset of 
asymmetry in massive progenitors of PNe. Further studies of the 
most obscured post-AGB sources must be pursued in the future. 
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